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Top flavour-changing neutral decays are extremely suppressed within the 
Standard Model (SM) by the GIM mechanism, but can reach observable rates in 
some of its extensions. We compute the branching ratios for t — > 07 and t — > eg 
in minimal SM extensions where the addition of a vector-like up or down quark 
\ singlet breaks the unitarity of the 3x3 Cabibbo-Kobayashi-Maskawa (CKM) 

\ matrix. The maximum rates obtained indicate to what extent present experi- 

mental data allow 3x3 CKM unitarity to be broken in these models, and are 
too small to be observed in the near future. As a by-product, we reproduce the 
, calculation of these branching ratios in the SM, and with an improved set of pa- 

, rameters we obtain values one order of magnitude smaller than the ones usually 

Ph. quoted in the literature. We study the CP asymmetries between the decay rates 

Ph. of the top quark and antiquark, which can be much larger than in the SM, also 

^ . as a consequence of the partial breaking of 3 x 3 CKM unitarity. 
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^ . 1 Introduction 

The arrival of top factories, LHC and TESLA, will bring a tremendous improvement 
in our knowledge of top quark properties |Lj, ||. In particular, the large top samples 
produced will allow to perform precision studies of top rare decays. In this field, 
flavour-changing neutral (FCN) decays t — > cZ, t — > 07, t — > eg, deserve special 
attention. Within the Standard Model (SM) they are mediated at lowest order in 
perturbation theory by penguin diagrams with quarks of charge Q = —1/3 inside the 
loop. Due to the smallness of down-type quark masses compared to Mw, these decays 
are very suppressed by the GIM mechanism, in contrast with processes like b — > 57, with 
diagrams with a top quark in the loop. This extra suppression results in decay rates 
O(10 -10 ) or smaller ||. On the other hand, in several SM extensions the branching 
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ratios for FCN top decays can be orders of magnitude larger. For instance, in two 
Higgs doublet models Br(t -> cZ) ~ 10" 6 , Br(t -> 07) ~ lO" 7 , Br(i -> eg) ~ can 
be achieved Q . In supersymmetric models with i? parity conservation these branching 
ratios can reach Br(i -> cZ) ~ 10~ 6 , Br(i -> 07) ~ 10~ 6 , Br(t -> c#) ~ 10~ 5 f| |. 

Here we are interested in the possible enhancement of these rates in models with 
vector-like quark singlets. The addition of quark singlets to the SM particle content 
represents the simplest way to break the GIM mechanism consistently. In these models, 
the 3x3 Cabibbo-Kobayashi-Maskawa (CKM) matrix is not unitary, and thus FCN 
couplings to the Z boson appear at tree-level. FCN couplings between light quarks 
are experimentally constrained to be very small, but this is not the case for the top 
quark. Actually, top FCN vertices can mediate the decays t — > uZ and t — > cZ, giving 
observable rates in models with up-type singlets Q. The largest branching ratios 
allowed by present experimental data are Br(£ — > uZ) = 7.0 x 10 -4 , Br(£ — > cZ) = 
6.0 x 10~ 4 0, much smaller than present direct limits Br(£ — > uZ), Br(t — > cZ) < 0.08 
d but still observable at LHC @ and TESLA |]| 0, In this Letter 
we investigate the enhancement of the branching ratios for the two other FCN decays, 
t — > 07 and t — > eg, in the presence of either up or down singlets. We find the rates 
of these processes allowed by present experimental constraints, and study how the 
GIM suppression takes place in these models. For completeness we also quote without 
discussion the results for t — > and t ug, which in the SM are suppressed by the 
ratio \V u b/V c b\ 2 with respect to the former, but in these SM extensions can have the 
same magnitude. 



2 Overview of the Lagrangian 

A full discussion of the Lagrangian in the weak eigenstate and mass eigenstate bases can 
be found for instance in Refs. Jl6|, 17]. Here we only collect the terms of the Lagrangian 
in the mass eigenstate basis relevant for our study. We consider the SM extended with 
n u up singlets and rid down singlets, with n u , rid arbitrary for the moment. The charged 
current Lagrangian is 

C w = -JLu L YVd L W+ + h.c. , (1) 

with V the generalised CKM matrix, of dimension (3 + n u ) x (3 + rid). The neutral- 
current Lagrangian describing the interactions with the Z boson is 

Cz = ~^ ( u ^X u u L - d L ^X d d L - 2s 2 w J£ u ) Z, , (2) 
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where X u , X d are hermitian matrices of dimension (3 + n u ) x (3 + n u ) and (3+rid) x (3 + 
no), respectively. These matrices can be related to the CKM matrix by I" = V V\ 
X d = V* V. The interactions with the unphysical charged scalars ^ are given by 

£ = —j3— u (M U VP L - VM d P R )d<j> + + h.c. , (3) 
V 2,Mw 

with Ai u and Ai d the diagonal mass matrices for the up and down quarks. The terms 
corresponding to the unphysical neutral scalar x are 

C x = ^-[u(M u X u P L -X u M u P R )u 

-d [M d X d P L - X d M d P R ) d] x • (4) 
Finally, the terms describing the interactions with the Higgs boson are 
C " = 7^^(M u X u P L + X u M u P R )u 



+d (M d X A P L + X d M d P R ) d] H . (5) 



In our analysis of t — > cy and t —>■ eg we discuss the two simplest cases: n u — 1, 
rid — (which will be called Model I) and n u = 0, rid — 1 (Model E). These two cases 
correspond to CKM matrices of dimension 4x3 and 3x4, respectively, and in both 
models the CKM matrix is a submatrix of a unitary 4x4 matrix. 



3 Calculation of the decay rates 

Using unbroken SU(3) x U(l) gauge invariance and the facts that final state particles 
are on-shell and the photon has transverse polarisation, the transition amplitude for 
t — > C7 can be written with all generality as 

M(t^ n ) = u(p c )[za^q l/ (A 1 + B, l5 )]u(p t )e;(q), (6) 

with p t and p c the momenta of the top and charm quarks, respectively, q = pt — p c the 
photon momentum and e its polarisation vector. This expression also assumes that 
the top quark is on-shell, which is an excellent approximation. A similar structure is 
valid for t — > eg, with form factors A g and B g . In order to compute the amplitude 
the form factors are written in terms of Passarino-Veltman integrals |18|] using FORM 



1~9|. The integrals are numerically evaluated using LoopTools ||20|| . The Feynman 



diagrams relevant for t — > cy in the SM and Model E are depicted in Fig. [I]. in the 
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SM the down-type quarks di in the loops are di = d, s, b, while in Model E there is 
an extra heavy quark B. The contributions of these diagrams to A 1 and B 1 in the 
't Hooft-Feynman gauge are collected in the Appendix. The diagrams relevant for A g 
and B g are the analogous to (la) and (lb) in Fig. |l| but replacing the outgoing photon 
by a gluon. 



W 4> 




(2a) (2b) 




(2c) (2d) 

Figure 1: Feynman diagrams contributing to the t — > cy decay amplitude in the SM 
and Model E. 



In Model I there are extra diagrams with up-type quarks Ui = u, c, t, T in the loops 
(see Fig. |^). These diagrams have one FCN vertex for Ui = c,t and two for Ui = u, T, in 
which case they are very suppressed. The flavour- diagonal vertices are modified with 
respect to the SM value. For instance, the diagonal couplings of a quark q = Ui, di to 
the Z boson are 
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as can be seen from Eq. (|2[), with the plus (minus) sign for up (down) quarks. The inter- 
actions with the unphysical scalar x an d the Higgs boson can be read from Eqs. (§J^). 
The contributions of these diagrams to v4 7 and _B 7 can be found in the Appendix. 





(3a) 



(3b,3c) 



Figure 2: Additional Feynman diagrams contributing to the t — > cy decay amplitude 
in Model I. The diagrams for t — > eg are similar, replacing the photon by a gluon. 



We perform the computation keeping all quark masses. For external quarks we use 
the pole masses mt = 174.3 GeV, m c = 1.5 GeV. For internal quarks it is more adequate 
to use MS masses at a scale 0(m t ), rather than pole masses. This is an important 
difference due to the strong dependence on the b quark mass as a consequence of the 
GIM suppression. For a pole mass mj = 4.7 ± 0.3 GeV, m h {rn t ) = 2.74 ± 0.17 GeV 



In the limit m c = the vector and axial form factors are equal: 



B-f, Ag 



Since m r is small, A, 



Bj, Ag ~ 



Bg and the effective couplings are predominantly 



right-handed. One important feature is that the form factors acquire imaginary parts 
from the contributions with d, s, b quarks (and u, c quarks in the extra diagrams 
present in Model I). These imaginary parts are one of the ingredients needed in order 
to have CP asymmetries T(t — > cy) ^ T(t — > cj), which will be analysed in detail later. 

From Eq. (H), the partial widths of these processes are 
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(8) 



with Cf = 4/3 a colour factor. In the SM, as well as in our models, the total width 



is dominated by the leading decay mode t — > bW + , T(t — > bW + ) ■ 
m t = 174.3 GeV, My/ = 80.39 GeV. The branching ratios are then 

r(* -> cry) o u ^ ^ c 9) 



Br(t -> cry) 



T(t -> bW+) 



Br(t -> eg) 



T(t -> bW+) 



1.57\V th \ 2 for 



(9) 
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We do not use the next-to-leading order partial width F(t — > bW + ) = 1.42 \ V t b\ 2 for 
consistency, because our calculation for t — > 07, t — > eg is at leading order. 

We have checked that using the set of input parameters of Ref. our results agree 
with the results presented there. For the calculation of the branching ratios within the 
SM we take \V US \ = 0.2224, \V ub \ = 0.00362, \V cb \ = 0.0402. These values are obtained 
performing a fit to the six measured CKM matrix elements, using 3x3 unitarity. The 
phase 5 in the standard parameterisation []22] is 5 = 1.014, obtained with a fit to e, 
e'/e, a^K s and \5m B \ (see Ref. ||). The SM predictions are 

Br(t^c 7 ) = (4.6^±0.4±J!)xlO- 14 , 

Br(*-c<7) = (4.6 ^ ± 0.4 ^-i) x 10-12 . ( 10 ) 

The first uncertainty comes from the bottom mass, the second from CKM mixing angles 
and the third is estimated varying the renormalisation scale between Mz (plus sign) 
and 1.5 m t (minus sign). These figures are ten times smaller than the ones quoted 
in Ref. ||, where the pole mass is used for the internal b quark (m?, = 5 GeV is 
assumed). The uncertainty in the top mass does not affect these values, because the 
partial widths of t — > cy, t — >• eg are proportional to 77t|, and the partial width of 
t — > bW + is approximately given by 



,2 ^3 



T{t^bW + ) = ^HW 2 ™ 4 



64tt' 101 



mf m? 



(11) 



Hence, the leading dependence on m t cancels in the ratios and the uncertainty in m t 
hardly affects the numbers quoted in Eqs. (lilf) . The SM predictions for t — ► wy and 
t — > ug are 

Br (t^ M7 ) = (3.7^±2.1^i)xl0- 16 , 

Br(f-u<7) = (3.7^|±2.1^; 5 7 )xl0- 14 , (12) 

suppressed by a factor | V^^/^cbl 2 — 8 x 10 -3 with respect to top decays to a charm 
quark. The uncertainties have the same origin as in Eqs. (|i~0D. 



4 CKM unitarity and GIM suppression 

Let us discuss how the GIM mechanism suppresses these processes within the SM and 
how this suppression can be partially removed with the addition of vector-like singlets. 
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We only study t — > cy, the discussion of t — > eg is formally identical. In the SM and 
Model I the form factors for the jtc vertex can be decomposed as 

i 

B, = £AsK)A^ (13) 

i 

where i — 1, 2, 3 in the SM and i — 1 • • -4 in Model E, f^A^i), f^Bi^i) are functions 
of the internal quark mass and A* t = V C iV t * are CKM factors. We have dropped the 
bar over m^, which are understood as MS masses. Since / 7 ^(mj) ~ /^(m;), we only 
analyse A T The mass dependence of the real and imaginary parts of / 7 a(^i) is shown 
in Fig. |. 




m. m. 



Figure 3: Loop functions f~,A( m i) and /^(mj) for down-type internal quarks (notice 
the different scales). 



In order to estimate the branching ratio for t — > cy we make the approximation 
f-yA( m d) = f-yA( m s) = /-^(O) for the moment. Then, using 

Kt + Kt + Kt + Af t = , (14) 

as implied by the unbroken row unitarity of the 3x4 CKM matrix V, we have 

A, = [/ 7 .H-/ 7 a(0)]^ + [/ 7 aN-/ 7 ,(0)]A« 

= /; A (m 6 )A^ + Um B )\» . (15) 

Therefore, the decay amplitude is actually controlled by the shifted function /^ A (mj), 
plotted in Fig. |3] as well. The parameter Xf t measures the orthogonality of the c 
and t rows of the 3x3 CKM submatrix Vs X 3 (see Eq. (0)), i- e. the breaking of 
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the GIM mechanism in this process. The SM limit is recovered setting, Af = 0, so 
the only contribution to the form factor is given by the small function fL A {m b ) ~ 
—9.1 x 1CT 9 — 4.7 x 1CT 9 i multiplied by X b ct ~ 0.04. With an extra down quark, there 
is a new term with a larger function fL A {rriB) — 4.9 x 10~ 7 + 6.0 x 10~ 6 z (for mg = 200 
GeV), which is however suppressed by Af . 

The parameter Af can be related to the breaking of the column unitarity of V^. 
This is easily understood, because if the columns of this submatrix are orthogonal, 
so must be the rows. The explicit relation can be written using the extension of the 

4 



Wolfenstein parameterisation [£3] in Ref. [|4j]. Assuming that Xd s , Xdb, X sb ~ A 

-bb 



1 — X„ ~ A 3 and 1 — X hh ~ A 3 , we have 



- Af t = E = X sb ~ ^X db + AX 2 (X bb - X ss ) - ^-X sb + 0(A 7 ) . (16) 
i=i 1 

This equation shows how the breaking of the orthogonality of the first three columns 

of V "propagates" to the second and third rows. The effect of the new quark can 

be estimated with A ~ 0.22, A ~ 1 and the typical values X db ~ 10~ 3 , X sb ~ 10~ 3 , 

X bb — X ss ~ 10 -3 H|, obtaining Af ~ X sb ~ 10 -3 . With this value the B contribution 

is 20 times larger than the b term, giving Br(t — > 07) ~ 10 -11 . 

In Model I, neglecting for the moment diagrams with two FCN vertices, A 1 can be 
decomposed as 

3 

^ 7 = fjA (m ) A* t + g lA X ct , (17) 
i=i 

with g lA the sum of the c and t diagram contributions, which is roughly of the same 
size as the f lA functions. In this model we have the relation 

Kt + Kt + Kt = X ct (18) 

expressing the non-orthogonality of the second and third rows of the CKM matrix, of 
dimension 4 x 3 in this case (compare with Eq. (|i"4l) ). Hence, the form factor is written 
as 

A 1 = [f lA {m b )-f lA {Q)}X b ct + [g lA + f lA {0)}X ct 

= f^ A (m b )X b ct + g' jA X ct , (19) 

with g' A = —4.4 x 10~ 6 — 4.8 x 10 -6 %. In this model the FCN coupling X ct can be 
X ct ~ 0.04 for V tb ~ 0.6 ||, yielding a branching ratio Br(t — > 07) ~ 5 x 10~ 8 . 

We note that the larger branching ratio achieved in Model I is not a consequence of 
the presence of a tree-level coupling Ztc, which appears in the expressions of the form 
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factors on the same footing as the parameter A^. Moreover, the loop integrals of the 



new physics contributions of Model I, g'^ A 

-6 , 



Vf 
-4.4 x 10' 



4.8 x 10" 6 i and of Model 



E, f' lA (mB) — 4.9 x 1CT 7 + 6.0 x 10 _ti i, are very similar. The only reason for the larger 
branching ratio in Model I is that X ct ^> A^, that is, unitarity of V3 X 3 can be broken to 
a lesser extent in Model E due to the strong requirements on FCN couplings between 
light quarks. Additionally, V t b can be smaller in Model I, and the total top width is 
reduced. 



We are also interested in the CP asymmetry 

r(* _> c 7 ) - r(t 



CTf) 



(20) 



1 Y(t -> cry) + T(t -> cry) 

This interest is mainly academic, because if the branching ratios are unobservable, even 
less are the asymmetries. However, the latter show how large CP asymmetries at high 
energy are possible in these SM extensions. Here we analyse in detail a 7 in Model E, 
the results for Model I are similar but more involved. The form factors for t — > cry are 



.4, 



B* 



and an analogous expression for B 1 . The asymmetry can be written as a 
with 



(21) 

N y /D 7 , 



-21m 



f',A( m b)f'U( m B) + f',B{ m b)f'*iB{ m B) 



Im 



B* 



|A^| 2 + 



+2 Re [/^K)/'*aK) + f^B( m b)f'* lB ( m B) 
A few comments are in order: 



\fU m B)\ 2 + \KB(m B )\ 
Re 



I A 



B\2 
ct\ 



\ b \ B* 

A ct A ct 



(22) 



1. The CP asymmetry a 7 is proportional to the imaginary part of the rephasing- 
invariant quartet A^ t A^* = V c bV* B V t * b VtB ■ This is expected from general grounds. 
In fact, it can be shown that in a model with an extra down singlet all CP 
violating observables at high energy (that is, when m Ui d,s ~ compared to the 
scale of energy involved) must be proportional to lmV c bV* B V t lV t B, lmVcbV* B XbB, 
Im VtbV tB X b B, or a combination of them |25| . 



2. The SM limit is recovered setting A^ = 0, obtaining a vanishing CP asymmetry. 



It is well known that CP asymmetries at high energy are very small in the SM p6 
due to: (i) the smallness of m u , and m s , what leads to a more efficient GIM 
cancellation; (ii) the small mixing between the top and the first two generations. 
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When all the quark masses are kept in the computation a non-vanishing but 
a negligible asymmetry a 7 ~ —6 x 10~ 6 is obtained. For the gluon case, the 
asymmetry a g ~ —5 x 1CT 6 is also extremely small. 

Since CPT invariance requires that the total width of the top and the antitop are 
equal, the different partial widths T(t — * cry) ^ T(i — > cry) must be compensated 
in other channel. In this case, the compensating decay channels are the SM 



leading modes t — > bW, t — > sW, t — > dW and their conjugate processes |26| . 



4. The large phases in the functions fy A , f' lB allow to obtain relatively large CP 
asymmetries, provided Im A^A^* is sizeable. 



5 Results 

We explore the parameter space of Models I and II to find the maximum values of 
Br(t — > try) and Br(t — > eg) allowed by present experimental measurements. The 
constraints on these models come from precision electroweak data, K and B physics 
and atomic parity violation (the details of the analysis can be found in Ref. ||). We 
take all the quark masses into account, and require that the mass of the new quark is 
larger than 200 GeV to satisfy the limits from direct searches. In Model I, assuming 
that the new quark has a mass rrtT = 200 GeV, we find the maximum rates 

Br(t try) = 4.5 x 10" 8 , 

Br(t^c^) = 8.9 xlO" 7 , (23) 

corresponding to \X ct \ = 0.037, \V t b\ = 0.58. (The branching ratios scale with |X ct | 2 
approximately.) For larger m^, the allowed values of \X ct \ are smaller 0, and these 
branching ratios decrease. The CP asymmetries corresponding to the figures in Eqs. (|23| ) 
are negligible, 

a 7 = -0.0006, 

a g = -0.002, (24) 

because the rates are dominated by the X ct term. The asymmetries can have values in 
the range —0.5 < a 7 < 0.4, —0.9 < a g < 0.6, but only reach the boundaries of these 
intervals for branching ratios much smaller than those in Eqs. (|23|) . The results for 
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decays to up quarks are a little larger, 

Br(t -> «7) = 4.6 x 10~ 8 , 

Br(t^wg) = 9.2xl0~ 7 . (25) 
In Model E, assuming that the mass of the new quark is mg = 200 GeV, we have 

Br(t -> try) = 4.5 x 10~ 12 , 

Br(t -> eg) = 6.6 x 10~ n . (26) 

These numbers are almost insensitive to the mass of the new Q = —1/3 quark for 
m>B > 200 GeV, as can be seen from Fig. [I], and show a small increase with tub- The 
corresponding asymmetries are 

a 7 = —0.05 , 

a g = 0.56. (27) 

In general, the CP asymmetries take values in the intervals — 1 < a 7 < 1, — 1 < a 9 < 1. 
The decay rates to up quarks are larger, 

Br(t -> uj) = 7.4 x 10~ 12 , 

Br(t -> up) = 9.5 x lO" 11 . (28) 
The branching ratios in Eqs. (|23|-28D are too small to be measurable in the near 



future. The estimated 3a sensitivities of LHC to these decays are Br(t — > cry) 
1.2 x 10" 5 0, Br(t -> «7) = 3.0 x 10~ 6 @, Br(t -> eg) = 2.7 x 10" 5 and Br(i 



ug) = 4.1 x 10 -6 ||28|| , with an integrated luminosity of 100 fb _1 . The TESLA sensitivity 
to t —>■ cry is better but not enough, Br(£ — > cry) = 3.6 x 10 -6 with a centre of mass 
energy of 800 GeV and a luminosity of 500 fb _1 . Hence, we observe that in models with 
up-type singlets the rates for t — > qZ, q = u, c can be observable but not the rates for 
t — > Ql and t — > qg, which are four and three orders of magnitude smaller, respectively. 
This fact contrasts with the results for two Higgs doublet models or supersymmetric 
extensions of the SM, where the branching ratios for t — > cZ and t — > C7 are similar, 
and the branching ratio for t — > eg is one order of magnitude larger. This difference 
would allow for a consistency check of the models, should a positive signal of top FCN 
decays be discovered. 
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A Form factors for t cry and t —> eg 

The contributions to and _B 7 of diagrams (la) and (lb) with an internal quark di 
are 

Qig 2 e 1 

Ay.ia = l - — Y^V 2i V^{(m t + m c )Co + (2m t + m c )C 1 + (m t + 2m c )C 2 

+m t C u + (m t + m c ) C 12 + m c C 22 } , (29) 
Qig 2 e 1 

B y la = l - — — — rV^* {(m t - m c ) C + (2m t - m c ) G x + (m t - 2m c ) C 2 

2 lD7T z 

+m t C n + (m t - m c ) C 12 - m c C 22 ] , (30) 
Qi9 2 e 1 
4M^ 16tt 2 



+ (m t m 2 dt + m c m 2 ) C u + [m t (m 2 c + m 2 di ) + m c (m 2 + m 2 ,.)] C 12 



B i,ib = - 4 ^ 2 e i67r 2 ^ V 3* {-m c (m 2 - m 2 d .) C x + m t (m 2 c - m\) Q 



+ (m t m 2 c + m c m 2 d J C 22 j , (31) 
Qi9 2 _e_ 1 

-2 
W 

+ (m t m 2 di - m c m 2 ) C n + [m t (m 2 c + m 2 di ) - m c (m 2 + m|)] C 12 

+ (m t m 2 c - m c m\) C 22 j . (32) 

The C's are functions of the external and internal masses, C(m 2 ,0,ml, M^,m d .,Wi d .) 
in the notation of Ref. p0|| . For t — > eg the contributions to the form factors A g and 



B g can be obtained replacing e by g s and setting Qi — 1 in Eqs. (|29H32|). The terms 
from diagrams (2a)-(2d) are 

# 2 e 1 

-4 7 ,2a = ——r^y2iV^ i {(2m t + m c )C + (Am t + m c )C 1 + (m t -m c )C 2 

4 1071^ 

+2m t C n + 2{m t - m c ) C l2 ] , (33) 
g 2 e 1 

^ 7 ,2a = —-^V 2i V^ i {(2m t -m c )C + (4:m t -m c )C l + (m t + m c )C 2 

+2mtCu + 2{m t + m c ) C 12 } , (34) 
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i 7 ,2b 



g 2 e 



4M 2 , 16tt 2 



l —V 2i V^ {m c (m 2 - m 2 d J d + (m t m 2 d . + m c m 2 t ) C u 



+[m t (m 2 d . - m 2 c ) + m c (m 2 - m 2 ,.)] C 12 ] , 



B, 



7,2b - 



;V 2i V£ {-m c {m 2 - m 2 d .) d + (m t m| - m c m 2 ) C n 



i 7 ,2c - 



-^7,2c 
^7,2d 



ff 2 e 1 
4M^ 16tt 2 

+[m t (m| - m 2 ) - m c (m 2 - m 2 ,.)] C i2 } , 
V^Vg* {m c C + m c Ci + m c C 2 } , 

V 2i V 3i {-m c C - m c Ci - m c C 2 } 
V 2i V*{-m t C 2 } , 
V 2i V*{-m t C 2 } . 



9 2 e 1 
4 16tt 2 

9 2 e 1 
4 16tt 2 

4 16tt 2 
S 2 e 1 



(35) 

(36) 
(37) 
(38) 
(39) 
(40) 



4 167T 2 

Here the arguments of the functions are C(m 2 ,m 2 ,0,M^,m|,M^). In Model I the 
contributions from diagrams (3a) and (3b) with internal quarks Ui are 

Qi9 2 e 1 



i 7 ,3a 



4c 2 , 16tt 2 



{ c ll [-( m t + m c ) C - (2m t + m c ) Ci - (m t + 2m c ) C 2 



B, 



7,3a 



i 7 ,3b - 



-m t C u - (m t + m c ) C 12 - m c C 22 ] + (d LR + c RL )2m Ui [Co + d + Cs]} , 

^Ac 2 6 167r 2 {° %LL ' _< ^ m< ~~ m ^ C ° ~ ^ 2m< _ m ^ Cl ~~ ~~ 2mc " > ° 2 
-m t C u - (m t - m c ) C 12 + m c C 22 ] + (c % LR - c RL ) 2m Ui [C + C x + C 2 ]j , 

Qi9 2 e 1 
8M^ 16tt 2 



B 



7,3b - 



X«X« {-m c (m 2 - m 2 J d + m t (m 2 Ui - m 2 c ) C 2 

— {m t m 2 u . + m c m 2 t ) Cn - \m t (m 2 c + m 2 Ui ) + m c (m 2 + m 2 .) C 
-(m t m 2 c +m c m 2 u JC 22 } , 

{+m c (m 2 - mlj Ci + m^m^ - m 2 ) C 2 



12 



(41) 



Q^ 2 e 1 
8M^ 16tt 2 

-(m t m 2 Ut - m c m 2 ) Cn - m t (m 2 + m 2 .) - m c (m 2 + m 2 ,) 
-(m t m 2 c - m c m 2 Ui ) C 22 } • 



C 



12 



(42) 



The constants c l LL , c l LR and c^ L are products of left-handed and right-handed couplings 
between (c, i) and (i,t), 



= [X^-S~-si 



4 



Xft — Sit- I 



13 



--LR 



^RL 



-6 A -s 2 



6 V 
°it Q °w 



(43) 



with Sij — 1 if i = j and zero otherwise. For these diagrams the arguments of the 
functions are C(m 2 , 0, m 2 , M§, m^, 1^). Finally, the Higgs contribution in diagram 
(3c) reads 

Qi9 2 e 1 



7,3c 



8M^ 16tt 2 



X«Xf 3 { [2mt< + m c (m 2 + m 2 ; 



Ci 



+ 
+ 



m t (m 2 c + m 2 .) + 2m c m 2 . C 2 + {m t m 2 Ui + m c m 2 t ) Cu 

m t (m 2 + m 2 .) + m c (m 2 + m 2 J Ci 2 + (m t m 2 + m c m 2 J C22} , (44) 



5 



7,3c 



Q^ 2 e 1 
' 8M^ 16tt 2 



XlXl { [2m t m 2 Ul ~ m c (m 2 t + m 2 j] C x 



+ 
+ 



m t (m 2 + m 2 J - 2m c m 2 J C 2 + (m t m 2 Ui - m c m 2 ) C n 

m 4 (m 2 + m 2 J - m c (m 2 + m 2 j] C12 + {m t m 2 c - m c m 2 u J C 22 } , (45) 



where the functions are C(m 2 , 0, m 2 , Mj^, m 2 ^ m 2 J. We take M# = 115 GeV. For 
£ — > eg the extra contributions in Model I can be obtained replacing e by g s and 
setting Qi — 1. 
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